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ABSTRACT 

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected 
galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected 
galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain 
robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS 
galaxies, which have reliable optical counterparts with r < 19.4 mag and well-determined 
spectroscopic redshifts. The overlap region between the two surveys is 12.6 deg 2 ; the matched 
sample has a median redshift of z k, 0.2. The cross-correlation of GAMA and H-ATLAS 
galaxies within this region can be fitted by a power law, with correlation length r « 4.63 ± 
0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies 
within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA 
galaxies of &h/^g ~ 0.6. Combined with clustering measurements from previous optical 
studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos 
with masses comparable to that of the Milky Way. The correlation function appears to depend 
on the 250 /im luminosity, L250, with bright (median luminosity ^L 2 5o ~ 1-6 x 10 10 L Q ) 
objects being somewhat more strongly clustered than faint (^L 2 5o ~ 4.0 x 1O 9 L ) objects. 
This implies that galaxies with higher dust-obscured star formation rates are hosted by more 
massive halos. 



1 INTRODUCTION 

It is well known that V galaxies are the largest contributors to 
the present-day stellar mass density (e.g. Li & White 2009). It is, 
however, not clear how star formation is distributed across galax- 
ies and halos of different masses. Previous studies show that in the 
local universe star formation takes place preferentially in low den- 
sity environments (e.g. Lewis et al. 2002; Heinis et al. 2009). The 



most commonly used estimators of the star formation rate (SFR) 
are based on the UV continuum or Ha, H/3 or [Oil] emission lines 
(e.g. Brinchmann et al. 2004; Salim et al. 2007). These are all sub- 
ject to uncertain dust extinction corrections, and so can greatly un- 
derestimate the SFRs in dust-obscured regions. Mid- and far-IR 
observations, which are sensitive to the energy re-emitted by dust 
heated by young stars, are therefore an essential complement to UV 
and optical tracers of star formation. Such dust is heated to temper- 
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atures of around 20 — 40 K, emitting thermal radiation which peaks 
at wavelengths around 100 /im. IRAS measured the far-IR emission 
from bright galaxies but more recent surveys of dust emission have 
focussed on either mid-IR (ISO, Spitzer) or sub-mm (e.g. SCUBA) 
wavelengths, missing the peak in the dust emission, and therefore 
requiring uncertain extrapolations to infer total IR luminosities and 
hence dust-obscured SFRs. The launch of Herschel (Pilbratt et al. 
2010) has now opened up the study of the universe at far-IR wave- 
lengths (60 — 700 /im), spanning the peak of the dust emission 
from star-forming galaxies, and allowing robust measurements of 
the dust-obscured SFR. The Herschel- ATLAS (H-ATLAS) survey 
(Eales et al. 2010) will provide far-IR imaging and photometry cov- 
ering the wavelength range from 110 /im to 500 /im, over an area 
550 deg 2 , much larger than previous surveys at these wavelengths 
such as BLAST (Devlin et al. 2009). 

Analysis of clustering statistics provides a simple but powerful 
way to investigate environmental effects, in this case of the SFR of 
galaxies. In this paper we perform a preliminary clustering analysis 
of a 4x4 deg 2 field observed during the H-ATLAS science demon- 
stration phase (SDP). Previous analyses of the H-ATLAS (Mad- 
dox et al. 2010) and HerMES (Cooray et al. 2010) surveys have 
focused on angular auto-correlations, with no significant signal in 
the former case and a significant detection in the latter. Here, we 
consider spatial cross-correlations of far-IR and optical galaxies, 
which can be used to derive the clustering bias and hence the char- 
acteristic mass of the host halos. We analyze a sample of ~1000 
H-ATLAS galaxies which have reliable counterparts brighter than 
r < 19.4 mag in the SDSS and spectroscopic redshifts measured 
by the Galaxy and Mass Assembly (GAMA) survey 1 . Currently, 
the overlap region between the H-ATLAS and GAMA surveys is 
12.6 deg 2 , and the spectroscopic redshift completeness is 99.7% 
for galaxies with r < 19.4 mag. 

A full analysis of the spatial auto-correlation function of H- 
ATLAS galaxies is given in van Kampen et al. (2010). Here we 
instead measure the cross-correlation function of H-ATLAS and 
GAMA galaxies, a statistic that provides a more robust and ac- 
curate estimate of the clustering bias of the H-ATLAS galaxies. 
There are at least two reasons why this is so. Firstly, the sample 
of H-ATLAS in the relatively small SDP survey area is small. By 
contrast, the number of GAMA galaxies in this area exceeds that of 
H-ATLAS galaxies by a factor of ~ 10. Secondly, the redshift dis- 
tribution of the GAMA galaxies can be robustly measured from the 
full GAMA survey (rather than from just the restricted SDP area) 
and for the estimator we employ knowledge of the H-ATLAS red- 
shift distribution is not required. Thus the systematic uncertainties 
due to cosmic variance are reduced. As a result, the estimate of the 
cross-correlation function of the relatively sparse H-ATLAS sam- 
ple with the more populous GAMA sample has much better statis- 
tics than the estimate of the H-ATLAS auto-correlation function 
alone. Finally, even though our sample is relatively small, using 
the cross-correlation technique allows us to investigate the depen- 
dence of clustering on far-IR luminosity by dividing the H-ATLAS 
sample into two subsets according to 250 /im luminosity. In this 

1 GAMA will eventually provide a highly complete, wide-area spectro- 
scopic survey of over 400,000 galaxies with sub-arcsecond optical/near-IR 
imaging (from SDSS, UKIDSS, VST, VISTA), and complementary obser- 
vations from the UV (GALEX) through to the mid and far-IR (WISE, HER- 
SCHEL) and the radio (ASKAP, GMRT). GAMA has so far surveyed 144 
deg 2 and the catalogue contains 95,000 galaxy redshifts to r-band magni- 
tude 19.4 with a redshift completeness of 98.7% (Driver et al. 2010, 2009; 
Baldry et al. 2010; Robotham et al. 2010; Hill et al. 2010). 



manner, we determine the clustering bias and infer the typical halo 
mass for each subset. 

Throughout this paper we assume a flat ACDM cosmology 
with fi m = 0.25, n A = 0.75 and H ( , = 73 km s" 1 Mpc" 1 . 



2 SAMPLE SELECTION 

We use data obtained by the Spectral and Photometric Imaging 
Receiver (SPIRE, Griffin et al. 2010; Pascale et al. 2010) in the 
16 deg 2 H-ATLAS science demonstration field 2 . In total there are 
6, 878 sources over an area 14.4 deg 2 that are brighter than the 
5a detection limit in one or more of the 3 SPIRE bands: 250 /im, 
350 /im and 500 fim (Rigby et al. 2010). The corresponding flux 
limits are 33, 36 and 45 mJy/beam. Below we work with the 
250 /im flux limited sample as this is the most sensitive band, has 
the best positional accuracy and was used for source detection in 
the catalogue that was matched to GAMA (Smith et al. 2010). 

A significant fraction of these 6,878 Herschel galaxies lie at 
low redshifts and have optical counterparts in the SDSS imaging 
catalogue. Sources with S/N>5 at 250 /im (6,621) were matched 
to the r-band selected (r <22.4) SDSS catalogue by Smith et al. 
(2010) using a likelihood ratio analysis (Sutherland & Saunders 
1992; Ciliegi et al. 2003) with a maximum 10" search radius. This 
leads to 4,756 sources which have at least one candidate optical 
counterpart in SDSS. A reliability value (7?lr) is then assigned 
to each of the optical candidates, which quantifies the probability 
that the counterpart is a genuine match. We discard candidates with 
Rlr. < 0.8 to remove unreliable matches, leaving 2,424 reliably 
matched sources. The angular overlap of the GAMA 9-hour field 
with H-ATLAS is not perfect and this reduces the survey region 
(hereafter GAMA-SDP) from 14.4 deg 2 to 12.6 deg 2 . Within this 
region there are 2,143 reliably matched sources. The spectroscopic 
redshift coverage of GAMA in this region is complete at 99.7% 
for an r-band Petrosian magnitude (corrected for Galactic extinc- 
tion) brighter than 19.4 mag. Imposing this cut leaves 969 galaxies 
which have measured spectroscopic redshifts and form the sample 
we analyse below (the H-ATLAS sample). A statistical analysis of 
the excess number of close pairs shows that 16% of GAMA sources 
brighter than 19.4 mag have a Herschel-ATLAS counterpart, of 
these, ~80% have directly identified reliable matches (Smith et al. 
2010). We only have spectroscopic redshifts for H-ATLAS galax- 
ies that have reliable matches in the r-band limited GAMA survey. 
Hence while we believe we a have complete representative sample 
of the local H-ATLAS galaxies we could in principle be missing 
galaxies which are bright at 250 /im but too faint for detection in 
the r band. This possibility can not be ruled out until we have spec- 
troscopic redshifts selected in the sub-mm. 

To k-correct the observed Herschel fluxes to the rest-frame 
250 pim, we assume that the dust emission has an SED of the form 

L„ ocB„(T)v , (1) 

where B V (T) is the Planck function. There are two parameters in 
this formula, the dust temperature, T, and the emissivity index, /3. 
We adopt the values, T = 28 K and (i = 1.5, derived by Amblard 
et al. (2010) by fitting to nearby H-ATLAS galaxies detected in at 
least three far-IR bands with a significance greater than 3a. 

The luminosity distribution at 250 [im (L250) is shown in 



2 PACS data, Ibar et al. 2010, are also available but are not used here. 
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Figure 1. Observed distribution of Herschel 250 (im rest-frame luminosities 
for the 969 galaxies well matched to SDSS galaxies with r < 19.4 mag. The 
dotted line corresponds to the threshold (2.5 X 10 24 W Hz -1 ) used to split 
the sample into bright and faint subsets (see text Section 2). 



Fig. 1. It peaks at around L 25 o = 3.2 x 10 24 W Hz , cor- 
responding to the local L* galaxies found by Dye et al. (2010). 
We further split the H-ATLAS sample into two subsets (indi- 
cated by the vertical dotted line in Fig. 1): bright sources with 
L250 > 2.5 x 10 24 W Hz -1 (corresponding to total IR luminosity, 
Lir = 5.0 x 1O 1O L0, based on Eqn 1, integrating from 8 -1000 
/im), and faint sources with L250 < 2.5 x 10 24 W Hz -1 . The faint 
subset consists of 484 galaxies and the bright one of 485 galaxies. 
The median values of L250 for the faint and bright H-ATLAS sub- 
samples are 1.3 x 10 24 and 5.0 x 10 24 W Hz -1 respectively (corre- 
sponding to total IR luminosities of 2.5 x 10 10 and 7.9 x 1O 1O L ), 
so that they differ by a factor 3 in typical luminosity. Fig. 2 shows 
the number counts as a function of the 250 fim flux. Although these 
two subsets are well distinguished in luminosity, they have similar 
distributions of observed 250 /im flux. 

The separation of the two samples by L250 is somewhat 
blurred by the uncertainties in the flux measurements and assumed 
k-corrections. Perturbing the luminosities according to the flux 
measurement errors in the H-ATLAS catalogue (Rigby et al. 2010) 
makes little difference with just 5% of the sample switching from 
the bright to the faint subsets. The k-correction depends on the val- 
ues of T and assumed in equation (1). The sample of Amblard 
et al. (2010) spans the ranges T=28±8 K and /3=1.4±0.1. This un- 
certainty can also scramble the luminosity subsets somewhat but 
even the most extreme choice of T = 36 K and /3 = 1.5 only 
switches 8% of the sample from the bright to the faint subsets. We 
return to the effect this might have on our clustering results in Sec- 
tion. 3. 

The distributions of apparent and absolute r-band magni- 
tudes (corrected for Galactic dust extinction) are shown in Figs. 3 
and 4 respectively. The r-band absolute magnitudes have been k- 
corrected to z = (Blanton et al. 2003). The r-band absolute mag- 
nitude for the full H-ATLAS sample peaks around —21.7, some- 
what brighter than the Milky Way. For comparison, we also in- 
clude the corresponding properties of the full GAMA sample in the 




log(S/mJy) 

Figure 2. Flux distribution of Herschel sources at 250 fim. Blue and red 
dashed curves are for the bright and faint Herschel sources respectively, 
with the black line showing the combination of the two. 



same sky area in Figs. 3 and 4. It can be seen that while there are 
more GAMA than H-ATLAS galaxies, their distributions of appar- 
ent and absolute r-band magnitude are similar. The extra galaxies 
in the GAMA catalogue may correspond to early-type and some 
late-type galaxies, for which the current SFRs are very low, leading 
to their absence from the far-IR survey. More detailed work on the 
properties of these galaxies is needed in the future. 

The redshift distributions of our samples are shown in Fig. 5 
as histograms. In each case, the upper histograms and curves cor- 
respond to the GAMA sample and the lower ones to the H-ATLAS 
sample. As expected, the more luminous H-ATLAS galaxies tend 
to lie at higher redshifts. The redshift distributions of the lumi- 
nous and faint galaxies cross at z ~ 0.2, which is roughly the 
median value for all the 969 H-ATLAS sources. To help interpret 
the cross-correlation of the faint and bright H-ATLAS sources with 
GAMA galaxies, we want subsets of the GAMA galaxies with sim- 
ilar redshift distributions to the corresponding H-ATLAS samples. 
To achieve this, we split the GAMA sample at M r = —21.2 mag 
into faint and bright subsets. It can be seen in Fig. 5 that this choice 
of dividing magnitude results in the corresponding subsets of H- 
ATLAS and GAMA samples having very similar redshift distribu- 
tions. The full, faint and bright GAMA samples have median abso- 
lute magnitudes M r of —21.5, —20.5 and —22.0 mag respectively. 



3 CORRELATION FUNCTIONS 

In this section, we first calculate the auto-correlation functions of 
the GAMA and H-ATLAS galaxies, then their cross-correlation, 
and finally the clustering bias of the H-ATLAS galaxies. The auto- 
correlation of the GAMA galaxies is needed for calculating the 
relative bias from the cross-correlation, while the H-ATLAS auto- 
correlation provides a consistency check on the results from the 
cross-correlation, and also allows us to compare with the auto- 
correlation results of van Kampen et al. (2010). 
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Figure 3. The r-band number counts per magnitude per square degree ver- 
sus apparent magnitude. As in Fig. 2, the blue and red dashed curves are for 
the blight and faint H- ATLAS sources, and black is the combination of the 
two. To compare, the green curve gives the number counts of all GAMA 
galaxies in the same region. 




-20 
M r [mog] 

Figure 4. Similar to Fig. 3, but as a function of r-band absolute magnitude. 
The curves are colour coded in the same way as in Fig. 3. 



3.1 Auto-correlation functions 

In this subsection, we estimate the auto-correlation function of the 
GAMA and H-ATLAS SDP samples and, for the H-ATLAS sam- 
ple, the dependence of clustering strength on the Herschel 250 
luminosity, L250- We begin by considering the correlation function 
in redshift space, £ (rx , ru), where r± and rii are the comoving sep- 
arations perpendicular and parallel to the line of sight respectively, 
and integrate this over the line-of-sight separation, rii, to obtain the 



Figure 5. Redshift distributions of each of our samples, plotted as the num- 
ber of galaxies per unit redshift per square degree. The colour coding is 
as in Fig. 3: red curves are for the H-ATLAS-faint and GAMA-faint sam- 
ples, blue curves are for the H-ATLAS -bright and GAMA-bright samples, 
black curves are for the full H-ATLAS sample and green for the full GAMA 
sample. For H-ATLAS, the histograms are the observed distributions, and 
smooth curves are fits to these distributions. For GAMA, the histograms 
show the redshift distribution of galaxies in the GAMA-SDP field, while 
smoothed curves are the fits to the full GAMA area. The fitted curves are 
collected in the bottom-right panel. 



projected correlation function. This removes the effect of peculiar 
velocities on the estimate of the spatial correlation function. 

There are several estimators for the auto-correlation function 
in the literature, all which require the generation of a uniform ran- 
dom catalogue with the same mask as the galaxy catalogue itself. 
In this work we adopt the estimator proposed by Hamilton (1993), 

DD(r ± ,r ]] )RR(r_ L ,r ]l ) 



[DR(r ± , n )]2 



(2) 



where DD(r±, ru), DR(r±,ru) and RR(r±,r^) are counts of 
data-data, data-random and random-random pairs, respectively. To 
generate smooth redshift distributions for the random samples we 
fit their redshift distributions with the functional form 



N(z) 



' exp(—flz v ). 



(3) 



The fits to the redshift distributions of GAMA and H-ATLAS 
(sub)samples are shown as smooth curves in Fig. 5. To obtain a 
robust estimate of the mean redshift distribution of the GAMA 
galaxies we made use of the full 144 deg 2 of the GAMA catalogue 
(~9xl0 4 galaxies), rather than just the subset that overlaps with 
the H-ATLAS area (~7 x 10 3 galaxies). The completeness mask of 
Norberg et al. (2010) was used to generate the random catalogue 
corresponding to the GAMA sample. 

Following standard practice, we estimate the projected corre- 
lation function, w(r p ), by integrating Eqn. 2 along line of sight ru : 



w(r p ) = w(r±) 



£(r_L,r||)dr||. 



(4) 



In reality, we cannot integrate to infinity. Instead, we have chosen to 
integrate to ±50 Mpc, but we test the impact of varying this limit. 
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Errors are estimated using the Jackknife technique. We split each 
galaxy sample into 16 equal area regions and then calculate the 
correlation functions for data taken from any 15 of these 16 regions. 
The scaled scatter of the Jackknife samples gives an estimate of 
the errors on the corresponding correlation functions (e.g. Norberg 
et al. 2009). 

The projected correlation function is related to the real-space 
correlation function by a simple Abel transform (Peebles 1980). 
For a power law, w(r p ) = Ar v ~^ ', the 3-D correlation function, 
£(r), is also a power law, £(r) = (r/ro) -7 . The parameters are 
related by 



r 7 



^T(7/2) 



r(i/2)r[( 7 -i)/2] : 



(5) 



where T(x) is the standard Gamma function. 

The two-point projected auto-correlation functions are plotted 
in Fig. 6 as red curves. To test the convergence of the line-of-sight 
integral in Eqn. 4, we show with green curves (here and later also 
in Fig. 7) the result of extending the integration out to 100 Mpc. 
The projected correlation function is seen to be insensitive to the 
precise choice of integration limit. 

For the GAMA-SDP sample (top row), the projected corre- 
lation functions are measured in the region of overlap with IT- 
ATLAS. The reason for remeasuring the GAM A correlation func- 
tions in this restricted area rather than showing the less noisy esti- 
mate from the full GAMA dataset (Norberg et al. 2010) is that we 
are interested in the relative clustering of H-ATLAS and GAMA 
galaxies and this choice will reduce the impact of sample vari- 
ance on the comparison. The GAMA-SDP correlation function can 
be well fitted with a power law (black dot-dashed curve). Fitting 
Eqn. 4 to the data in the range (1 - 12) Mpc and using Eqn. 5, we 
find r = 5.96 ± 0.62 Mpc and 7 = 1.87 ± 0.21. This is consistent 



with the values for L* optical galaxies (with M" 



-21.1) es- 



timated in the SDSS: r = 6.6 ± 0.3 Mpc and 7 = 1.87 ± 0.03 
(Zehavietal. 2010). 

The lower panels of Fig. 6 show the projected auto-correlation 
functions for our three H-ATLAS samples. The best-fit values of 
ro and 7 for these are summarised in Table 1. van Kampen et 
al. (2010) have carried out a more detailed analysis of the auto- 
correlation function of H-ATLAS galaxies using the angular corre- 
lations in redshift slices. They obtain a best estimate of the spatial 
clustering length, averaged over the redshift range 0.1 < 2 < 0.3, 
of ro = 5.5 ± 0.9 Mpc. Our estimate of ro given in Table 1 for 
the auto-correlation of our full H-ATLAS sample is consistent with 
this. 

For comparison, the best-fit power law for the GAMA-SDP 
auto-correlation function for the full GAMA sample is reproduced 
by a grey line in all of the other panels of Fig. 6. The full H-ATLAS 
sample is somewhat less clustered than the full GAMA-SDP sam- 
ple. The faint H-ATLAS galaxies appear to have similar clustering 
to the full H-ATLAS sample, while the bright H-ATLAS galaxies 
appear to be more strongly clustered. However, the statistical uncer- 
tainties in the estimates for these small samples are clearly rather 
large and, moreover, systematic errors could be introduced by fit- 
ting smooth curves to their noisy redshift distributions. These limi- 
tations are largely overcome in the next section where we measure 
the clustering of the H-ATLAS galaxies by cross-correlating with 
the much larger GAMA-SDP sample. Furthermore, by estimating 
the GAMA-SDP radial selection function using the full GAMA 
survey covering an area about 10 times larger than the GAMA- 
SDP region, systematic uncertainties in the modelling of the radial 
selection function are significantly reduced. 



In Fig 6 it is apparent that our jackknife error bars are some- 
times noisy as witnessed e.g. by the large error bars at~0.7Mpc 
or >6Mpc for the GAMA bright sample, or by the small error 
bars on the GAMA faint sample on scales below 2Mpc. Further 
investigation has revealed that this is a result of our small sample 
and occurs because the clustering on particular scales can be dom- 
inated by one or two structures and so vary significantly in just 
one or two of our jackknife samples. Such fluctuations are smaller 
for our cross-correlation samples, discussed below. Thus the er- 
rors quoted for the correlation length, ro, for the bright and faint 
GAMA auto-correlation samples have significant uncertainty, but 
the cross-correlation results and their error bars are more robust. 
The diagnostic tests used for the robustness of the clustering errors 
are similar to those presented in Norberg et al. (2010). 

3.2 Cross-correlation functions 

The cross-correlation function in redshift space of H-ATLAS with 
GAMA galaxies is estimated using 



£(rx,ni) = 



HG(r ± ,r n )RR(r ± , ril 
HR(r ± ,r n )GR(r ± ,r l{ 



1, 



(6) 



where HG, HR, GR, RR are counts of H-ATLAS-GAMA, H- 
ATLAS-random, GAMA-random and random-random pairs, re- 
spectively. In each case, the random sample is generated so as to 
match the redshift distribution of the GAMA galaxies. Thus, for 
our estimates of the cross-correlation functions, at no point do we 
need to fit the noisy redshift distributions of the small samples of H- 
ATLAS galaxies. As for the auto-correlation functions, we calcu- 
late the projected two-point cross-correlation functions according 
to Eqn. 4 and estimate the errors using the jackknife technique. 

The projected cross-correlation functions are shown in Fig. 7. 
The top panel shows the GAMA-H-ATLAS result when the limit 
of integration in Eqn. 4 is taken to be 50 Mpc (red curves) and 
100 Mpc (green curves) respectively. The dot-dashed line is the 
best fitting power law to the 50 Mpc estimate. It shows that the 
H-ATLAS-GAMA cross-correlation function is well fitted by a 
power law, with ro = 4.63 ± 0.51 Mpc and 7 = 2.05 ± 0.31, in- 
dicating that the clustering of the H-ATLAS galaxies is weaker 
than that of GAMA-SDP galaxies. This inferred difference between 
the strength of the H-ATLAS and GAMA-SDP clustering appears 
larger than suggested by comparing the upper and lower left panels 
of Fig. 6, or the values of ro in Table 1. This might be the result 
of a bias in the redshift distribution of the random samples for the 
H-ATLAS galaxies, which is obtained by fitting a smooth function 
to noisy data (§ 3.1). 

The lower panels in Fig. 7 show cross-correlation functions for 
subsets of luminous and faint H-ATLAS and GAMA galaxies. For 
comparison, this best-fit line to the GAMA-H-ATLAS function is 
replicated in grey in these panels. Again, we find that the clustering 
of faint H-ATLAS galaxies is weaker than that of the bright galax- 
ies. The estimates of ro and 7 for these samples are summarized in 
Table 1. 

As discussed in Section 2, there are uncertainties in the 
250 /im k-correction and the flux measurements. Adopting the 
most extreme perturbation to the k-corrections (T = 36 K and 
/3 = 1.5, see Section 2) and perturbing the fluxes according to 
the measurement errors quoted in Rigby et al. (2010) shifts the ro 
values of our estimated H-ATLAS autocorrelation functions by an 
amount comparable to the quoted 1-ct statistical uncertainty. This 
variation is largely caused by the limited size of these samples and 
the resulting uncertainty in fitting their redshift distrbutions. The 



© 0000 RAS, MNRAS 000, 000-000 



6 Qi Guo et al. 




Figure 6. Two-point projected auto-correlation functions. From left to right and top to bottom, they correspond to the GAMA, GAMA-faint and GAMA-bright; 
H-ATLAS, H-ATLAS-faint and H-ATLAS-bright samples respectively. The red curves show the result of truncating the line-of-sight integration in Eqn. 4 at 
50 Mpc and the green dashed curves at 100 Mpc. The black dash-dotted lines are power-law fits to the data in red for full GAMA and H-ATLAS samples only. 
To aid comparison, the fit to the GAMA auto-correlation function for the full sample (top left panel) is repeated as a grey curve in the other panels. Error bars 
are estimated using the Jackknife technique. 



cross-correlations on which we focus, and which do not depend on 
the redshift distributions of the H-ATLAS samples, are much less 
affected by the uncertainties in the k-corrections and flux measure- 
ments. In this case, the same perturbations affect the ro values, by 
no more than 15% of their quoted statistical error and so make a 
negligible contribution to the uncertainty in our results. 



3.3 Bias of H-ATLAS galaxies 

To interpret the meaning of the estimated large-scale cross- 
correlation functions, consider the simple linear bias model in 
which the auto- and cross-correlation functions of H-ATLAS and 
GAMA galaxies are related to the auto-correlation function, £ m , of 
the mass at redshift z — by 



£ H (r) = bu(z) 2 D 2 (z)U.v) 
£ G (r) = b G {z) 2 D 2 {z)i m {v) 
£ HG (r) = bn{z)b G {z)D 2 {z)£, m {v), 



(7) 
(8) 
(9) 



where the subscripts H and G denote H-ATLAS or GAMA respec- 
tively, D(z) is the linear growth factor of the perturbations in the 
mass, and £ m (r) is the auto-correlation function of the dark matter. 
In this case, the projected cross-correlation function that we have 
estimated is related to that of the mass at z = through 



where the average product of the bias and growth factors is given 
by 

, n 2, » Jn G (z)nn(z)b G (z)b H (z)D 2 (z){^)dz 

(b G b H D [z)) = ,(11) 

J n G {z)nn{z) (^) dz 

and tih(z) and n G (z) are the mean space densities of the H- 
ATLAS and GAMA samples at redshift z. For the auto-correlation 
function of the GAMA galaxies this reduces to 



w G (r p ) = (b G D 2 (z))w m (r p ), 
where 

a a fn 2 G (z)b 2 G (z)D 2 (z) (^)dz 
(b 2 G D 2 (z)} = - J GV ' GV ' V ' KdzJ 



jn G (z)(^)dz 



(12) 



(13) 



wn G (r p ) 



{bub G D (jz))w m (r p ), 



(10) 



The relative bias of the H-ATLAS and GAMA galaxies is then, 

&hg = w Gil (r p )/w GG (r p ) = (bnb G D 2 (z)) / (b 2 G D 2 (z)) . (14) 

In principle, this depends on both the bias parameters &h, b G and 
on D(z). However, since by construction the redshift distributions 
of the full/faint/bright H-ATLAS samples match well with those 
of the corresponding (full/faint/bright) GAMA samples, the depen- 
dence on D(z) will approximately cancel. If the bias parameters 
&h and &g evolve with redshift in the same way, then this evolution 
will also approximately cancel out in the relative bias. This is the 
reason why we cross correlate H-ATLAS faint/bright with GAMA 
faint/bright instead of all GAMA galaxies. 
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We estimate the mean relative bias of H-ATLAS and 
GAMA galaxies using, 



7rel S ^HG 
&HG = — 



El/a? 



(15) 



u 

Q. 



where 6^q' is obtained directly from the measured projected H- 
ATLAS-GAMA cross-correlation function and the GAMA auto- 
correlation function (rather than from the fits given in Table 1), and 
o"i represents the Jackknife error on 5^' estimated at each pair 
separation. This simple estimator ignores correlations between the 
measurements at different separations and so may not be optimal, 
but we do take account of such correlations in estimating the error 
on 6hg- O ur error on 6hg ls estimated using the Jackknife tech- 
nique, by calculating the mean bnGj f° r eacn Jackknife sample , 
(assuming the same values of a\ as used in Eqn. 15) and then look- '""q, 
ing at the scatter in values between Jackknife samples. 

Our estimates of &5Jq are shown in Fig. 8. For the full H- 
ATLAS and GAMA samples, the mean relative bias over the range 
of separations 2-8 Mpc, where the two-halo term dominates and 
where we have good statistics, is b^ G (alV) = 0.61 ± 0.08. Thus, 
we conclude that the clustering strength of H-ATLAS galaxies is 
significantly weaker than that of GAMA SDP galaxies. This impor- 
tant conclusion is revealed only by taking advantage of the cross- 
correlation function technique. As shown in Table 1, our estimates 
of the auto-correlation functions are much too noisy (and proba- 
bly subject to systematic errors) to detect any difference between 
the two galaxy samples. From the cross-correlation of the faint H- 
ATLAS with the faint GAMA samples, we obtain a relative bias 
of faint) — 0.67 ± 0.13, while from the cross-correlation 
of the bright H-ATLAS with the bright GAMA galaxies, we obtain 
bnaibright) = 1.04 ± 0.22. 

To convert the estimates of relative bias into values of the ab- 
solute bias for the different H-ATLAS samples, we need to know 
the absolute bias of the different GAMA samples. For this we use 
the results of Zehavi et al. (2010), who measured the clustering as a 
function of r-band luminosity in the SDSS, and combined that with 
a theoretical prediction for the clustering of the dark matter in the 
ACDM cosmology. An important qualification is that the values of 
bias measured by Zehavi et al. effectively apply at the average red- 
shift of the SDSS, z ~ 0.1. The bias of r-band selected galaxies is 
expected to evolve with redshift, but quantifying the size of this ef- 
fect for the redshift range z < 0.5 probed in the present paper must 
await a detailed clustering analysis of the full GAMA redshift sur- 
vey. Here, we will simply assume that the bias factors for GAMA 
and H-ATLAS galaxies can be taken to be constant over the redshift 
range studied here. We therefore use Eqn. (10) from Zehavi et al., 
scaled to as = 0.8, to calculate the value of the bias as a function 
of r-band absolute magnitude. 

Our full, faint and bright GAMA samples have median abso- 
lute magnitudes M° A = -21.3, -20.3 and -21.8 respectively (k- 
corrected to z — 0.1 to be consistent with Zehavi et al.), implying 
average r-band bias factors of be, = 1.17, 1.05 and 1.29. This then 
leads to absolute bias values of 6 H = 0.71 ± 0.09, 0.70 ± 0.14 and 
1.34 ± 0.28 respectively for the full, faint and bright H-ATLAS 
subsamples. We find that the bright H-ATLAS galaxies are more 
strongly clustered than the H-ATLAS population as a whole at the 
2a level, which confirms the trend seen from the H-ATLAS auto- 
correlation functions in Fig. 6. This result implies that the excess 
clustering of the bright H-ATLAS galaxies reflects a genuine and 
strong dependence of clustering on far-infrared luminosity and thus 
on star formation rate. We detect no significant difference between 
the bias of the faint H-ATLAS galaxies and that of the population as 
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Figure 7. Two-point projected cross-correlation of all H-ATLAS with all 
GAMA galaxies (top), faint H-ATLAS with faint GAMA galaxies (middle) 
and bright H-ATLAS with bright GAMA galaxies (bottom). As in Fig. 6, 
red curves show the result of intergrating Eqn. 4 to 50 Mpc and the green 
dashed curves to 100 Mpc. Black dash-dotted curves show the power-law 
fits to the red curves over the range 1-12 Mpc. For comparison, the fit to 
the all H-ATLAS-GAMA cross-correlation function is replicated as grey 
lines in the lower two panels. Error bars are estimated using the Jackknife 
technique. 



a whole. This result, however, could be affected by our assumption 
of a constant bias over the redshift of interest. 

The final step is to use the estimated clustering bias of H- 
ATLAS galaxies to constrain the masses of the halos hosting them. 
In the ACDM model at the present day, the halo bias is a very weak 
function of halo mass for halos less massive than 10 12 M Q and in- 
creases rapidly with increasing halo mass at higher masses (Mo & 
White 2002). Using the fitting formula for bias as a function of halo 
mass at z = from Seljak & Warren (2004), obtained from sim- 
ulations of a ACDM universe, we infer an average host halo mass 
log 10 M/M Q » 12.1 i^ 5 (or a 2a upper limit log 10 M/M Q < 
12.8) for the full H-ATLAS sample. We find very similar values 
for the faint H-ATLAS subsample, log 10 M/M Q » 12. 0l^ 71 (or 
a 2a upper limit log 10 M/M Q < 13.0). For the bright H-ATLAS 
sample, the average halo mass is log 10 M/M@ « 13.6lg'4. The 
more luminous H-ATLAS galaxies thus appear to be hosted in sig- 
nificantly more massive halos than the faint ones. Note that given 
the large errors in the estimates of halo masses, it is reasonable that 
the 2a upper limits on the host masses of the faint and bright sub- 
samples are both higher than that of the full sample. 
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Table 1. The correlation length, ro, and slope 7, for the power-law fits to our auto- and cross-correlation functions and the mean redshift, z mcan , number of 
galaxies, -/V ga i , and relative bias of each sample. 
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Figure 8. The relative bias foj^Q estimated from the ratio of the projected H- 
ATLAS-GAMA cross-correlation function, u>HG( r p). to the corresponding 
GAMA auto-correlation function wq (r p ) . The black curve is for the bias 
of the full H-ATLAS sample relative to the full GAMA sample, while the 
blue curve is for the bright H-ATLAS galaxies relative to bright GAMA 
galaxies, and the red curve is for faint H-ATLAS galaxies relative to faint 
GAMA galaxies. Errors are estimated using the Jackknife technique. 



4 CONCLUSIONS 

We have used a subset of the H-ATLAS galaxies in the SDP field, 
which have spectroscopic redshifts from the optical GAMA red- 
shift survey, to calculate the projected cross-correlation functions 
of far-IR and optically selected galaxies. We find that these H- 
ATLAS galaxies (which have a median redshift z ~ 0.2, median 
250 fim luminosity L250 ~ 2.5 x 10 24 W Hz -1 , and median total 
IR luminosity Lir ~ 5.0 x 10 10 Lo) are significantly less strongly 
clustered than the optically selected GAMA galaxies (which have 
a median absolute magnitude, M r = —21.5 mag) at the same red- 
shifts. This effect is also seen (though with lower significance) in 
the auto-correlations of the H-ATLAS and GAMA galaxies. 

From the cross-correlation analysis, combined with the pre- 
viously measured clustering of optical galaxies in the SDSS, we 
find that H-ATLAS galaxies are less clustered than the dark mat- 
ter, with an average bias b — 0.71 ± 0.09. This implies a typical 
host halo mass of ~ 1.25 x 10 12 M Q for the H-ATLAS galaxies in 



our sample (which are mostly at low redshift), comparable to the 
halo of the Milky Way. These preliminary results for the host halo 
masses of the H-ATLAS galaxies are consistent with the theoretical 
predictions of Lacey et al. (2010), who find a typical halo mass of 
1.6 x 10 12 M Q . (Note that Lacey et al. (2010) used the halo bias for- 
mula of Sheth et al. (2001) which predicts a somewhat larger bias 
than the Seljak & Warren (2004) formula used here at low masses.) 

We also split our H-ATLAS sample into subsamples of high 
and low far-IR luminosity, and investigate their clustering proper- 
ties. Both the cross- and auto-correlation functions suggest a de- 
pendence of clustering on far-IR luminosity over the range Lir = 
2.5 x 10 10 — 7.9 x 1O 1O L , with the bright galaxies being more 
strongly clustered than the faint ones at 2a significance, imply- 
ing that the more luminous galaxies are hosted by more mas- 
sive dark halos. The average halo mass for the bright sample is 
around 4 x 10 13 Mq and the 2<r upper limit for the halos hosting 
the faint sample is 1O 13 M0. The dependence of clustering on far- 
IR luminosity that we find here appears significantly stronger than 
the model predictions of Lacey et al. (2010) who find A/haio ~ 
1.3 x 10 12 and 2.0 x 1O 12 M for galaxies of comparable lumi- 
nosities to our faint and bright subsamples. It will be interesting 
to test whether this discrepancy persists in the full H-ATLAS sur- 
vey. As luminosity and redshift are correlated in a flux limited sam- 
ple, our high L250 luminosity subset has a higher median redshift 
than its fainter counterpart. Hence, in principle, strong evolution of 
clustering with redshift could be contributing to our inferred de- 
pendence of clustering on luminosity. We will be able to directly 
address this ambiguity with the much larger full H-ATLAS sam- 
ple by splitting the sample into redshift bins. When completed, this 
survey will enable comprehensive investigations of the clustering 
and environments of star-forming galaxies. 
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